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ABSTRACT 

Physical constraints of any real system can have a drastic effect on its performance. Some of the more recognized 
constraints are actuator and sensor saturation and bandwidth, power consumption, sampling rate (sensor and control- 
loop) and computation limits. These constraints can degrade system’s performance, such as settling time, overshoot, 
rising time, and stability margins. In order to address these issues, researchers have investigated the use of robust 
and nonlinear controllers that can incorporate uncertainty and constraints into a controller design. For instance, 
uncertainties can be addressed in the synthesis model used in such algorithms as or // . There is a significant 
amount of literature addressing this type of problem. However, there is one constraint that has not often been 
considered; that is, actuator authority resolution. In this work, thruster resolution and controller schemes to 
compensate for this effect are investigated for position and attitude control of a Low Earth Orbit formation flight 
system. 

In many academic problems, actuators are assumed to have infinite resolution. In real system applications, such as 
formation flight systems, the system actuators will not have infinite resolution. High-precision formation flying 
requires the relative position and the relative attitude to be controlled on the order of millimeters and arc-seconds, 
respectively. Therefore, the minimum force resolution is a significant concern in this application. Without the 
sufficient actuator resolution, the system may be unable to attain the required pointing and position precision control. 
Furthermore, fuel may be wasted due to high-frequency chattering phenomena when attempting to provide a fine 
control with inadequate actuators. 

To address this issue, a Sliding Mode Controller is developed along with the boundary Layer Control to provide the 
best control resolution constraints. A Genetic algorithm is used to optimize the controller parameters according to 
the states error and fuel consumption criterion. The tradeoffs and effects of the minimum force limitation on 
performance are studied and compared to the case without the limitation. Furthermore, two methods are proposed to 
reduce chattering and improve precision. 


I. INTRODUCTION 

Physical constraints can be enforced in hardware or software. Software limitations ensure the life of the system by 
reducing wear and tear on the actuators. For example, a very large step size or impulse control command might 
damage the equipment, and the saturation of the actuator may degrade the system performance and stability 8 . The 
effects and control of other hard nonlinearities, such as saturation and bandwidth limits have been investigated by 
numerous researchers. For instance, References 5-7 discuss how limitation on the maximum force and slew rate can 
be incorporated into the synthesis model of a or /i controller. Nonlinear controllers such as sliding mode can 
also be used successfully with actuator and slew rate saturations as discussed in Refs. 2-4. In sliding mode control, 
the saturation is addressed by the control command saturation and the sliding surface together. The slew rate is 
accounted for by the controller’s parameter. 

A physical constraint that is not routinely considered is actuator resolution and its effect on the performance. In low- 
precision applications, actuator resolution is not a concern since the desired precision may be orders of magnitude 
above the actuator capability. In a high-precision formation flight application, thruster resolution (minimum thrust 
impulse) must be addressed. Lack of actuator resolution can lead to performance degradation, such as high 
frequency chattering. For example, controller designs based on the sliding mode concept are susceptible to high 
frequency chattering unless an integrated boundary layer function is implemented. For cases in which the thruster 
has a minimum force limitation, high frequency chattering may still exist even with the boundary layer integrated 
into the sliding mode controller design. This is also addressed in this paper. 

Section II lists the symbols, abbreviations and acronyms used in this paper. Section III discusses the dynamics 
model, the control strategy, and the optimization method on which the actuator limitation effect and proposed 
control methods are based. This work is implemented in simulation on a Low Earth Orbit formation mission. 
Section IV discusses the effects of the thruster limitation in terms of predefined performance criteria, fuel 
consumption, control precision and chattering phenomenon. Sections V and VI present two methods that were 
studied in order to mitigate chattering and improve control precision in the presence of the thrust resolution issue. 
Section VII provides conclusions and discussions for future work. 
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II. SYMBOLS, ABBREVIATIONS AND ACRONYMS 


() : the time derivative with respect to a localized 
frame 

()’ : the true anomaly derivative 
| || : norm 2 
| | : absolute value 

mN: milliNewton 
LEO: Low Earth Orbit 
BLC: boundary layer control 
diag : diagonal matrix 
GA: genetic algorithm 
SMC: sliding mode control 
/ : true anomaly 


J : inertia matrix of the satellite 
p : relative position vector 
R : absolute position vector 
R : magnitude of the absolute position 
u : normalized force w.r.t. the mass 
p : gravitational coefficient 
X U ) : j th derivative of * 

LVLH: local vertical local horizontal coordinate 
frame 

Subscript F or L: follower or leader satellite 
Subscript p: position 

Superscript L: expressed in the leader’s local frame 


III. DYNAMICS, CONTROL, OPTIMIZATION 2 ’ 3 ’ 4 , AND SIMULATION CONDITION 


Dynamics Model 

This paper adopts the generalized coupled dynamics model 2 ,3, 4 for the leader-follower satellites pair in the 
formation flight system as the basic model. Eq. (1.1) shows the generalized relative position model expressed in 
LVLH and can be used in both the circular and elliptic reference orbit. 



Figure 1 LVLH Coordinate of the Leader’s Frame 
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The detailed information of the dynamics can be found in Ref [2] and is not presented in this paper. 


( 1 . 1 ) 


Sliding Mode Controller 

If x = [/?, , p 2 , p,] r is the system state, then the dynamics model of the leader-follower pair can be expressed as 


( 1 . 2 ) 


where B = is the input influence matrix and / = describes the dynamics of the system. The sliding 

surface is defined as 


dt 


(1.3) 
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Here x = x-x DESIRED is the tracking error, and = ^desired ~~ • Then the SMC is given by Eq. (1.4), in 

which, x ( r n ~ l) = [jft* xj.^ _I) ••• and ^sgn^)^,- sgn(s,)] mxl . k is defined as the sum of the uncertainty in 

the model plus j\. If the model is assumed to be exact, then k { = 77 , . Based on the stability analysis from the literature, 
A,. > 0 and 77. > 0 9 . In sliding mode control, parameter k t affects x andx , whereas, \ affects X . 

u = B-\x { r l) -f-ksga(s)] (1.4) 

In order to mitigate the chatter, the boundary layer control (BLC) concept is introduced as in Refs. [1, 3, 4]. In the 
BLC, the sliding surface includes a hard switching outside of the boundary, and a soft switching inside of the 
boundary. With the addition of BLC, the performance of the controller is smoother and more efficient. Equation 
(1.5) is used to replace the sign function in Eq. (1.4), where s i is a nonzero positive value. When e t goes to infinity, 
the switching becomes the hard switching. Smaller £■ produce a smoother but slower response. 

sgnfs,) s t > 1 

$ 

signiSj)^ 5flt(— ) = < s, , / = 1,2,.. .777 (1.5) 

S. - s,< 1 

l 


Parametric Optimization 

Nine parameters are used to design the SMC-BLC controller for the relative position model (Eq. ( 1 . 6 )). For a 
system that is highly nonlinear and has a long evaluation time, the conventional large-scale optimization does not 
work well. Therefore, a genetic algorithm search approach is utilized in order to expedite the selection of the 
parameters and also to obtain an “optimized” set of parameters. The population size and number of generation are all 
twenty. The bit length is forty, and the crossover and mutation probabilities are 0.9 and 0.05 separately. In order to 
reduce the calculation time, thruster working frequency is chosen as 0.5 Hz, and each evaluation time equals the 
period of the leader’s orbit. Since the full parameter space is not considered, there are no claims of global optimality. 

In SMC: k eR 3xl A eR 3xl 

-+p -*p 

In BLC: s p eR 3xl ( 1 . 6 ) 

IV. EFFECTS OF THRUSTER LIMITATIONS 
Simulation Condition 

The LEO formation system considered has both the leader and follower in 500 km altitude circular orbits. The 
leader’s orbit is equatorial and the follower’s orbit has an inclination of 0.1°. At simulation epoch, the leader’s true 
longitude is 0° and the follower is -0.1°. The small angular difference ensures the relative position of the formation 
system is on the order of 10 km and that the relative angles are small enough for small angle approximation. It is 
assumed that the simulated system has a thruster configuration that provides uncoupled thrust and torques in all 
directions. The objectives are maintaining the satellites in the fixed formation and not exceeding the physical 
limitation. The design requirements are to produce a stable system and find a balance between control precision, 
chatter, and fuel consumption. In future works, specific requirements such as jitter, controller bandwidth, attitude 
and position accuracy, and slew requirements will be used. From Section IV through Section VI, the relative 
position performance will be used to show the effects of the thruster limitation, and comparison of the two proposed 
chatter mitigation methods. 

Thruster without Limitation 

In order to generalize this work, the command and thruster forces are normalized with respect to the mass of the 
spacecraft. This also allows for force equivalence between different size spacecrafts in the formation. In the first the 

case, which does not use the maximum and minimum limits, control parameters k = [100,1 00,520] r , A = [0.1, 0.1,0. if , 

s = [0.01,0.01, 0 .Olf are chosen by trial and error. Figure 2 and Figure 3 illustrate the force commands where 
maximum and minimum thruster limitations do not exist. Because x and y directions are in plane and z is out of 
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plane, the initial velocity and force command in the z direction is much larger than that in the x and y directions. The 
saturation level in the z direction is due to the saturation in the sliding surface instead of the saturation in the thruster. 



Tim© Seconds Tim© Seconds 


Figure 2 Normalized force commands in X and Y Figure 3 Normalized force commands in Z direction 
directions of LVLH of LVLH 

Thruster with Maximum Limitation 

The thruster configuration is assumed to provide different maximum forces in each direction [1 mN 2 mN 2 N] in 
xyz respectively. Figures 4 through 6 display the simulation results when the thruster has the maximum limitations. 
If thruster doesn’t have a minimum thrust limitation, the resulting relative position error is on the order of 
micrometer. Figure 7 provides the results of the relative position error in z direction when no maximum limitation is 
enforced. A comparison of Figure 5 and Figure 7 illustrates that the saturation in the thruster induces a longer 
transient stage, more oscillations, and larger overshoot. 



Tim© Seconds Tim© Seconds 

Figure 4 Relative position in x of LVLH with Figure 5 Relative position in y of LVLH with 

maximum force limitation maximum force limitation 



Time Seconds Tim© Seconds 

Figure 6 Relative position in z of LVLH with Figure 7 Relative position in z of LVLH without 

maximum force limitation force limitation 
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Thruster with Maximum and Minimum Limitation 

Simulations are conducted and compared for different maximum and minimum force ratios. The maximum force 
ratio is with respect to F_5=[lmN 2mN 2N]. For example, a maximum force ratio of two corresponds to 2 *F_s 
means the maximum force is 2*[lmN 2mN 2N]. The minimum ratio is with respect to the maximum force so 2% 
means that the minimum force the thruster can provide is 2% of the maximum force. The performance criterion is 

calculated by y = 1| x(t) - x 4es (r)|ft , where x Des represents the desired state vector. The limitations constraints is 

defined as 


0 if\u\<=u m 

Signer™ '/|«|>=“m 
u other 


(1.7) 


Table 1 demonstrates that the minimum thrust ratio does not significantly affect the performance because the 
performance criterion is dominated by the transient stage and the control precision is in the stable stage. Instead, the 
performance criterion is reduced as the maximum thruster force increases. The increase in thrust levels reduces the 
settling time. 


The fuel consumption of the different cases is compared Table 2. For the relative position and velocity, error is 
linearly correlated with the control effort in SMC + BLC, therefore, the trend in Table 2 is similar to that of Table 1. 
Considering the maximum force level, performance criterion, and fuel consumption in both tables, the maximum 
ratio of 4 is chosen as the base line thruster in later simulations. Therefore, the maximum forces the thruster can 
provide in three directions are [4 8 8000] milliNewton. 


Table 1 Performance Criterion Comparison for Different Maximum and Minimum Level 


Min Ratio -» 

0% 

2% 

4% 

6% 

8% 

10% 

Max Ratio 4 

1 

15.1416 

15.1416 

15.1417 

15.1421 

15.1426 

15.1405 

2 

4.4977 

4.4979 

4.4988 

4.4994 

4.5005 

4.5012 

4 

1.4726 

1.4749 

1.4812 

1.4900 

1.5005 

1.5109 

6 

0.8161 

0.8250 

0.8459 

0.8735 

0.8900 

0.9221 

8 

0.6872 

0.7043 

0.7381 

0.7947 

0.8479 

0.8194 

10 

0.6872 

0.7128 

| 0.7689 

0.8226 

0.8194 

0.9820 


Table 2 Fuel Consumption in sum sense / = f 7 ( V \u^dt 


i=x t y,z 


Min Ratio -> 

0% 

2% 

4% 

6% 

8% 

10% 

Max Ratio 4 

1 

0.1083 

0.1082 

0.1081 

0.1081 

0.1081 

0.1079 

2 

0.0622 

0.0620 

0.0619 

0.0618 

0.0617 

0.0616 

4 

0.0397 

0.0395 

0.0393 

0.0391 

0.0389 

0.0388 

6 

0.0335 

0.0332 

0.0329 

0.0327 

0.0324 

0.0321 

8 

0.0316 

0.0312 

0.0308 

0.0304 

0.0300 

0.0296 

10 

0.0316 

0.0310 

0.0304 

0.0300 

0.0296 

0.0292 


Based on the tables above, if the limitation constraint technique is used for the same thruster, then the performance 
criterion and fuel consumption do not change significantly. In Figures 8 - 15, the time responses of the relative 
position error and control efforts are shown. The precision of the relative position is increased with the thruster 
resolution capabilities (i.e., better precision is obtained with the 0% case), as seen in Figures 8-11. In the z 

direction, the control precisions are on the order of 0(10~ n ) , 0(1(T 3 ) , and 0(1 (T 3 ) km, respectively for the 0%, 4%, 
and 10% cases. It should also be noted that higher control precision is achieved at the expense of significant 
chattering. 
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Figure 8 Relative position error in x of LVLH for 


different min force level 



Figure 10 Relative position error in z of LVLH for 
different min force level 



Figure 9 Relative position error in y of LVLH for 
different min force level 



Time Seconds 

Figure 1 1 Relative position error in z of LVLH for 
different min force level from 3000 to 6000 seconds 


Figures 12 - 15 illustrate the time response of the thruster forces. The 0% (infinite resolution) thruster case is 
smooth and has no chattering, whereas the 10% case gives more chatter. Figure 15 contains the time response for the 
10% case, which has several spikes due to the resolution. In each Figure there are 3 minimum force percentages. 
The higher the percentage, the more likely chatter will occur. 




Figure 12 Normalized thruster force in x of LVLH Figure 13 Normalized thruster force in x of LVLH 

from 0- 400 seconds 
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Figure 14 Normalized thruster force in y of LVLH Figure 15 Normalized thruster force in z of LVLH 


V. COMPARISON OF DIFFERENT PREFILTERS 

As seen in the results from the previous section, the chattering and lowest achievable relative position error are 
associated with the discontinuity of the thruster force commands. In this section, we discuss a pre-filter concept to 
mitigate this discontinuity generated by the SMC. The idea is to attempt to map the control commands generated by 
die SMC to “smooth” thrust commands. A total of three pre-filters are compared. Pre-filter 1 is the baseline case 
used normally and is described in Eq. (1.8); pre-filters 2 and 3 are described in Eqs. (1.9) and (1.10), respectively. 
Figure 16 provides a graphical description of all three pre-filters. It can be observed that pre-filters 2 and 3 have 3 
discontinuities and pre-filter 1 has 4 discontinuities. Both pre-filters 2 and 3 try to balance the total forces by 
increasing the force below the minimum, and decreasing the force above the minimum. 


Pre-filter 1 (cut off case): 


0 

sign(u)u m 


'/ M >= “»x 

other 


Pre-filter 2: 


s/gw(w)[i^ + 


u 


“max X rati0 

sign(u)u m 


-(«» 


ifu = 0 

-O] if «*0|w| <=«„,„ 
if M>«m» 


( 1 . 8 ) 


(1.9) 


Pre-filter 3: 




Sfe»(«)[- 


'-u +u„ 


sign(u)u a 


if u = 0 

] if u*0and\u\<=u mMX 

if 


( 1 . 10 ) 


Figures 17-22 present the relative position error and resultant force for the three pre-filters. Although the pre-filter 
1 has a larger stable-state error, which is on the order of centimeter level, it results in very low chattering. In Table 3, 
for the fuel consumption case, pre-filter 1 also yielded the lowest fuel consumption, whereas, the other performance 
criteria for all three cases were similar. Pre-filters 2 and 3 both generate the force above the required one, therefore, 
chattering is more severe than the pre-filter 1. 
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Prefilter 



Figure 16 Three pre-filters 



Figure 18 Relative position in Y of LVLH 



Figure 20 Normalized force in X of LVLH 



Figure 22 Normalized force in Z of LVLH 



0.25 


N» 
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■g 0.05 1 
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Figure 17 Relative position in X of LVLH 

Comparison of Different Prefitter Effects 


- 0.1 



Prefilter 1 

Prefilter 2 with ratio 1.1 . 

Prefilter 3 

k 

■ 

r . 


500 1000 1500 2000 2500 5000 


Time Seconds 

Figure 19 Relative position in Z of LVLH 



Figure 21 Normalized force in Y of LVLH 
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Table 3 Performance and Fuel Consumption Comparison for Different Pre-filters 



Performance Criterion 

Fuel Consumption in 
Sum Sense 

Pre filter 1 

1.5109 

0.0388 

Pre filter 2 

1.4730 

0.2473 

Pre filter 3 

1.4800 

0.2456 


VI. COMPARISON OF DIFFERENT BLC 

The boundary layer function is modified in order to mitigate the chattering. In Figure 23, three boundary layer 
functions are introduced which have 2, 5, and 3 discontinuities respectively. 


BLC 1: 

sign(s) 

s 

€ { 

sgn(s) = - 

s 

S 


. e ‘ 





Figure 23 Three different BLCs 


Table 4 shows the comparisons of the different boundary layer control schemes based on cost, fuel consumption, 
control precision and chattering. Here, the chatter is defined as the number of jumps in the control effort during 
simulation time period from 4000 to 6000 seconds in all directions. Cases 3 and 5 use the GA to determine the 


optimal (p^ . The cost is calculated by J = ^j[e T Qe+uRu]dt . In order to generate relatively equal contributions 


from all three directions in the cost, the weights are chosen to be Q = diag[2000, 1000,1] and R = diag[\ 000, 100,1] 


since the z-direction’s relative position error and control efforts are very different from those of the x- and y- 
directions. Fuel consumptions for all the cases are on the same order. However, without GA, the BLC 2 and BLC 3 
cases yield high cost functions and low precisions. The last column in Table 4 provides a qualifier of the chatter 
from 4000 to 6000 seconds of the simulation for three directions. 


It is worth noting that when the parameters for BLC 2 and BLC 3 are obtained from the GA, both cases converge to 
BLC 1 and the optimal results are the same as that of BLC 1. Therefore, it can be concluded that BLC 1 provides the 
best solution among all three BLC designs. The time responses of the relative position error and control efforts are 
the same as the pre-filter 1 case shown in Figures 17-22. 
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Table 4 Comparison of different BLC 



min 

Cost 

Fuel 

Control Precision in Km 

Chattering / 4000~6000 sec 

BLC 1 

0 

1.1423 

0.0388 

© 

© 

5 

i 

[2; 0; 0] 

BLC 2 

0.2 

423.4600 

0.0370 

[10- 2 ; lO' 2 ; lO' 2 ] 

[1; 8; 0] 

BLC 2 & GA 

1.22e-5 

1.1423 

0.0388 

[ 1 0" s ; 10" s ; 10“ 3 ] 

[0; 0; 0] 

BLC 3 

0.2 

371.5605 

0.0367 

[ 1 0- 2 ; 10' 2 ; 10~ 2 ] 

[1;7; 0] 

BLC 3 & GA 

6.65e-5 

1.1398 

0.0387 

[ 1 0 5 ; 10' 5 ; 10~ 3 ] 

[1;2;0] 


VIL CONCLUSION 

The addition of actuator saturation can degrade the transient performance of the system. The other extreme, actuator 
resolution limitations, can result in chatter, low fuel efficiency and a reduction in the achievable precision. In this 
paper, the resolution of the thruster is considered in the control problem of the formation flying system. It has been 
shown that when the thruster resolution is considered, chatter is introduced into the system and the performance of 
the system is degraded. Two methods were investigated to mitigate this problem. The first was an application of a 
pre-filter to smooth the command sent to the actuators. The pre-filter is used to compensate the discontinuity effects 
due to the minimum thrust level. 

Based on the fuel consumption consideration, the pre-filter 1 is the preferred one. If the chattering and control 
precision are considered more, the pre-filter 2 case is preferred 

Based on chattering and control precision results, pre-filter 1 is the preferred one. For this approach, three different 
prefilter designs were considered. The second approach was to modify the BLC; three modifications were also 
considered. The simulations, which were based on GA, illustrated the convergence of the three BLC schemes and 
that showed BLC1 has the best performance when considering the cost, fuel consumption, control precision and 
chattering. 

There are currently efforts in this endeavor; for example, the colloidal thruster being developed by Busek can 
provide very small force, on the order of 1 micro-N. However, while colloidal thrusters may be able to achieve the 
types of actuator resolutions needed, they may not be able to deliver the levels of maximum thrust needed. 

Currently, the maximum force from these types of thrusters are on the order of 25 micro-N which may be too small 
to satisfy the mission requirements. Further study is needed to define the trades between actuator precision and 
mission accuracy. 
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